Dissipative periodic waves, solitons and breathers of the nonlinear Schrodinger 

equation with complex potentials 
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Exact solutions for the generalized nonlinear Schrodinger (NLS) equation with inhomogeneous 
complex linear and nonlinear potentials are found. We have found localized and periodic solutions for 
a wide class of localized and periodic modulations in the space of complex potentials and nonlinearity 
coefficients. Examples of stable and unstable solutions are given. We also demonstrated numerically 
the existence of stable dissipative breathers in the presence of an additional parabolic trap. 
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I. INTRODUCTION 

Dissipative wave phenomena in nonlinear media with 
complex parameters are attracting nowadays a great deal 
of attention. They appear naturally in optics of active 
media like, for example cavities, active fibers, etc. [l| 
and in quantum mechanics when inelastic interactions of 
particles with external forces are accounted 0. More 
recently, there has been a particular interest in the dy- 
namics of nonlinear waves in periodic complex potentials 
due to possible applications to matter waves in absorb- 
ing optical lattices [!, H| and to periodic modulations of a 
complex refractive index in nonlinear optics 0, @ . In this 
last case special attention was devoted to the so called 
VT potentials (i.e. complex potentials for which the 
nonlinear Schrodinger (NLS) equation becomes invariant 
under the parity and time-reversal symmetry) Q , whose 
remarkable linear properties are presently intensively in- 
vestigated [1, II]. In particular, it has been shown that in 
nonlinear media with PT-symmctric linear damping and 
amplifications, stationary localized and periodic states 
may exist ■ The stability problem of these nonlinear 
structures is still an open problem. Also, it is not investi- 
gated so far whether similar structures with real energies 
could also exist in a general complex potential, thanks 
to the balance between nonlinearity, dispersion, gain and 
dissipation. 

In the present paper we shall address these problems 
with a twofold aim. From one side we derive a set of ex- 
act soliton solutions of the nonlinear Schrodinger (NLS) 
equation with complex potentials and show that the non- 
linearity management technique of the gain-loss profile 
can be an effective mechanism for providing stability. 

We find exact solutions by adopting an inverse en- 
gineering approach, 



developed in [11| for the case of 
the conservative NLS equation with periodic coefficients. 
More specifically we assume a specific pattern for the so- 



lutions and determine a posteriori the potentials which 
can sustain such pattern as a solution. Using this ap- 
proach we determine a general class of potentials which 
may support dissipative periodic waves and solitons in 
the form of elliptic functions. We show that some of 
these solutions my be stable with respect to small per- 
turbations. 

From the other side, we show that the derived solutions 
can be used to explore interesting dynamical behaviors 
of localized modes of the complex NLS equation. To 
this regard we investigate stable breather solutions which 
originate from localized exact solutions when a parabolic 
trap is switched on as additional real potential. Using the 
strength of the parabolic trap as a parameter we show 
the occurrence of a bifurcation from an attractor center, 
corresponding to a stationary soliton, to a limit cycle, 
corresponding to a stable breather mode. 

The stability of stationary periodic solutions, single 
hump solitons, and breathers opens the possibility to ex- 
perimentally observe such waves both in nonlinear optics 
and in BEC with optical lattices in presence of dissipative 
effects. 

We remark that although the determination of the po- 
tentials from the solutions is opposite to what usually 
occurs in practical contexts where potentials are a priori 
fixed, there are physical situations in which an inverse ap- 
proach can be experimentally implemented. An example 
of this is the case of a Bose- Einstein condensates nonlin- 
ear optical lattices which is produced by means of spatial 
modulations of the interatomic scattering length. In the 
mean field approximation the system is described by the 
NLS equation with spatially modulated nonlinearity of 
the type considered in this paper. Since the modulation 
of the scattering length (nonlinearity) in a real experi- 
ment is controlled by external magnetic fields via a Fesh- 
bach resonance, in principle one can construct arbitrary 
potentials to support specific prepared states, implement- 
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ing thus an inverse engineering approach. Another situ- 
ation, where the potential can be adjusted in accordance 
with the profile designed a priori, is the nonlinear optics 
of waveguide arrays, where the thermo-optics effect is 
employed using heaters producing the temperature gra- 
dients properly distributed in space. This gives the hope 
that some of the exact solutions reported in this paper 
could be indeed observed in experimentally contexts. 

The paper is organized as follows. In section II we 
introduce the model equation and explain the inverse 
method used to determine the solutions. In Section III we 
apply the method to the case in which there are only lin- 
ear real and complex potentials. As an example we derive 
solutions in the potential which has the 'PT-symmetry. 
We show that some of these solutions can be stable, what 
is confirmed by both the linear stability analysis and by 
numerical time evolutions of the complex NLS equation. 
In section IV we do similar studies for the case of non- 
linear lattices. In section V we show an example of exact 
stable solutions in correspondence of more general po- 
tentials, while in section VI we use this solution to nu- 
merically show the existence of dissipative breathers in 
the complex NLS equation with additional real parabolic 
trap. Finally, in section VII the main results of the paper 
arc briefly summarized. 



II. THE MODEL AND THE METHOD 

We consider the generalized NLS equation with varying 
in space complex linear potential, Ui(x) = Vi(x)+iWi(x), 
and complex nonlinearity parameter U n i = V„i(x) + 
iW ni (x) (hereafter V =ReU and W =Imf7) 

#t = —ifixx + <r|V»l V + WW + U n i(xM 2 ^, (1) 

where a = ±1 and is introduced in an explicit form in 
order to facilitate transition to the case of the homoge- 
neous nonlinearity (just putting V n i(x),W n i(x) = 0). In 
the particular case of the VT symmetry the invariancc 
imposes restrictions on the linear and complex poten- 
tials: namely Vi 7 V n i must be even and Wi 7 W n i must be 
odd functions of the space variable, respectively. 

We will be interested in the solutions of the form 
tp = A(x) exp (i[6(x) — ut]) , where A(x) and 9{x) are 
real amplitude and inhomogeneous phase of the mode. 
Substituting this ansatz in Eq. (TTJ) we obtain the system 
of equations 

uA + X -A xx - ^v 2 - a A 3 - VA - V nl A 3 = 0, (2) 
A x v + -Av x -WiA-W n iA 3 = 0, (3) 

where v = 6 X . 

Let us now assume that linear and nonlinear potentials 
Viix) and V n i (x) are given, and pose the problem of de- 
signing dissipative terms Wi(x) and W n i(x) to obtain a 



given solution A(x) fTTJ] . We emphasize that this choice 
is only for illustration proposes and any pair of the four 
functions vi, n i(x) and Wi, n l{x) can be chosen as a priori 
given. 

As first step, we prove that this is indeed possible, 
provided that A(x) is a bounded function [l3| satisfying 
the condition 

\f(x)\ < const, where f(x) = A xx /A (4) 

for all x. Indeed, after dividing @ by A we obtain and 
explicit expression for the velocity 

v 2 (x) = 2(w - Vi) + f(x) - 2(a + V n i)A 2 . (5) 

As it is clear this equation always has solutions for suf- 
ficiently large frequency u). Next, substituting ([5]) in (O 
we find the expression for the dissipative potentials 

If we want to consider only nonsingular dissipative terms, 
then we have to impose one more constrain on the field 
A, which must be considered together with (|4]) 

A x v , . 

— — < const, (7) 

and assume v to be a smooth function of x. 

In what follows we limit ourselves only by this kind 
of dissipation. This allows us to indicate immediately 
several types of admissible solutions: 

• Type 1: A is bounded and has no zeros. These, for 
example are functions like dn(a;, k) etc. 

• Type 2: A is bounded, has no zeros in any finite do- 
main, but decays as |x| — > oo. These are solutions 
like 1/ cosh(a;), exp(— x 2 ), etc. 

For these types of solutions condition (|7|) is automati- 
cally satisfied if we restrict our consideration to bounded 
conservative potentials. 

Further, we relax the requirements for A allowing it to 
have zeros but we require that condition (JT)) holds. Then 
the condition (Q| leads to another type of solutions: 

• Type 3: A has zeros which coincide with those of 
A xx . These are solutions of the type ^sn(a;, k), etc. 

All solutions reported below belong to one of these three 
types. We would like to notice here that A must not 
necessarily be taken in the form of elliptic functions. We 
choose elliptical functions just as an example represent- 
ing periodic solutions with a given parameter (the elliptic 
modulus). Alternatively we could choose the field A in 
the form of, for instance, Gaussian function. Another re- 
mark we would like to make here is that if A is a solution 
of the homogeneous NLS equation, then the coefficients 
Vi and V n i are constants and one has to determine Wj 
from either eq. ^ or ^ for a chosen W n i- Notice, 
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that this procedure implies some freedom, in the sense 
that one can also think of determining W„i for chosen 
Wi. However the latter might require some additional 
conditions on the linear dissipative potential Wi because 
sn-like solutions have zeros. In general terms this case is 
mentioned in the paper as a solution of Type 3. 

In closing this section we remark that from the above 
analysis it follows that for all constructed potentials 
one can define a function f(x) such that the nonlin- 
ear solutions are also solutions of the linear problem 
A xx — f{x)A = 0. This observation leads to an immedi- 
ate interesting conclusion. Assume that V n i(x) < and 
choose f(x) = — 2(uj — V{). Then, for a periodic functions 
Vi, the solution A is nothing but a linear Block func- 
tion of the linear potential 2Vi(x). Indeed, in this case 
A xx - 2V t (x)A = -2uA and v(x) = y / 2\V n i\A, which 
means that condition (|Tj) is satisfied. 



III. CASE OF LINEAR COMPLEX 
POTENTIALS 

According the scheme described above one can design 
the dissipation of the system in a way to provide the any 
desired field configuration supported by the given linear 
and nonlinear conservative potentials. While the zoo of 
available nonlinear patterns is not limited, physically the 
most relevant cases are those corresponding to the sta- 
ble solutions. Below we present several examples of the 
stable field distribution, completing the analysis with a 
contrasting behavior of unstable patterns. 

Before going into detail we also notice that the diver- 
sity of the physical parameters at hand makes the analy- 
sis of particular cases cumbersome. It turns out, however 
that at least one of them can be scaled out. In particu- 
lar, the amplitude of the nonlinear conservative potential 
can be scaled out by the renormalization of the amplitude 
of the solution. Therefore in what follows this value is 
chosen to be fixed. 

We start the analysis of the particular cases with the 
natural test example where A(x) is chosen to be a so- 
lution of the standard NLS equation, i.e. having con- 
stant real linear and nonlinear potentials. We respec- 
tively set Vi =const and V n i = 0. By means of straight- 
forward algebra one obtains, from ([5]) and (|6]), the re- 
lation Wi(x) = —W n i(x)A 2 . As it is clear this relation 
includes the particular case Wi = W n i = 0, which pre- 
cisely corresponds to the integrable NLS equation. Note, 
however, that the obtained relation allows for a general 
class of equations, all having the chosen function A(x) 
as a solution and the ratio of the linear and nonlinear 
complex parts of the potential fixed to —A 2 (x). 

Now we turn to the case where Ui (x) ^ and U n i (x) is 
a real constant, which will be chosen to be ±1, i.e. to the 
case where the inhomogeneity is only linear and nonlinear 
gain or loss are absent. The first example is a cnoidal 
wave (solution of Type 3 according to the classification 
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FIG. 1: (color online) (a) The distributions of the real (solid 
black) and imaginary (dashed red) parts of the linear potential 
given by Eq.([9| for k = 1 and Ao = 0.92, (b) the distribution 
of the amplitude |^| (solid black, left vertical axis) and the 
phase gradient v (dashed red, right vertical axis) of field ij). 
(c),(d) show the same but for k = 0.85 and Ao = 0.75. All 
panels are for Vo = —0.2 and Wo = 0.44. 

introduced in the previous section) 

A = A Q cn(x, k), (8) 

embedded in the linear lattice potential Vi(x) = 
Vocn 2 (a;, k) with a constant focusing nonlinearity cr = — 1. 
Requiring lo = 1/2 — k 2 , the respective pattern can be 
created with only the help of the linear dissipative term 

Wi = -W sn{x, fc)dn(x, k), W = ^ A 2 - V Q - k 2 , 

as it follows from ^ . The hydrodynamic velocity of the 
obtained solution is given by v = V2W Q cn(x,k). The 
phase is: 9 = (v / 2Wo/fc)arccos(dn(x, x)). The obtained 
complex potential Ui(x) is "PT-invariant. 

As it follows from © the nonlinear pattern A{x) ex- 
ists only for the amplitudes above the threshold A^' = 
yVa + k 2 and grows with the intensity of the dissipative 
part. If k =/= then the solution is periodic. However if 
k = 1 then our solution is localized and has no zeros. So 
if k = 1 then the solution belongs to Type 2. The typical 
field distributions for solutions are shown in Fig. [1] 

Let's now consider the stability of the solutions. To 
do this we linearize Eq. ([I} around the examined soliton 
solution. The solution of the linearized equation can be 
sought in the form tp(t,x) = J2 m ( Pm(x)exp(X m t)- The 
spectrum A of the small perturbations has both a con- 
tinuous part, associated to non- localized cigenfunction 
and a discrete part associated to localized modes (p m . 
The existence of A with positive real part means that the 
corresponding eigenmode will grow exponentially in time 
and, thus, the examined solution will be unstable. In this 
paper we study the stability numerically substituting the 
spatial derivatives by their discrete analogs (we used 5- 
point approximation). Then the problem of finding the 
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FIG. 2: (color online) (a) The spectrum for k = 1, A = 1.04, 
Vo = — 1 and Wo = 2.2 for the case when potentials are given 
by @ and the background solution is given by Eq.©. The 
arrows show the direction of the motion of the eigenvalues 
when Vo,; an d Wo,l increase (b) The eigenvector of the unsta- 
ble mode. 



eigenvalues A governing the stability of the solution is 
reduced to the diagonlization of a matrix. 

Discussing the stability we should first notice that if 
Vo = Wq = then Eq. (JTJ is nothing but the stan- 
dard NLS equation. The spectrum of small perturba- 
tion on the background of the Schrodinger soliton has 
two pairs of degenerate zero eigenvalues corresponding 
to phase and translational symmetries. The introduction 
of a nonzero Ui(x) preserves the phase invariance but 
breaks the translational symmetry. In Fig. [5] we show 
the evolution of the eigenvalues in the complex plane as 
the linear potentials are increased with the ratio Vq/Wq 
kept constant. 

We see that by increasing the potentials away from 
zero, the splitting of the zero eigenvalue present at Vq = 
Wo = gives a pair of eigenvalues lying in the gap of 
the continuum which move in opposite directions along 
the imaginary axis without producing any instability of 
the solution. We have therefore that for sufficiently small 
U(x) the soliton solution is always stable. However, when 
the pair of the eigenvalues reach the border and collide 
with the continuum they generate a quartet of complex 
eigenvalues with two unstable modes generating an oscil- 
latory instability of the soliton. 

The development of the instability is illustrated in 
Fig. [3J Panel (a) shows the instability of a soliton with 
the parameters Vo — — 1 and Wq — 2.2. The instability 
results in an infinitely growing localized peak in the area 
where the gain is positive. The formation of this peak is 
clearly seen on the picture. Panel (b) shows the decay 
of a periodic structure because of the instability present 
at the parameters Vq = —0.2, Wo = 0.44 (let us mention 
that at this parameters localized solution is stable). One 
can see that the instability is of the same kind and results 
in the formation of growing peaks in the areas of positive 
gain. 

The solution ijHJ) is periodic and describes a pattern 
with spatially alternating currents, such that the average 
hydrodynamic velocity is zero (see Fig.HJ. It is not diffi- 
cult to present an example where a pattern corresponds 
to a nonzero currents. To this end we consider the same 
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FIG. 3: (color online) (a) The decay of an unstable solitary 
solution ©, A = 1.04, Vo = -1 Wo = 2.2, k = 1 (b) The 
decay of an unstable periodical solution, Ao = 0.7, Vo = —0.2, 
Wo = 0.44, k = 0.8. 
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FIG. 4: (color online) (a) The distributions of the real (solid 
black) and imaginary (dashed red) parts of the linear potential 
given by Eq.© for k = 1, A = 0.71, (b) the distribution of 
the amplitude of \ip\ (solid black, left vertical axis) and the 
phase gradient (dashed red, right vertical axis) of field ip given 
by Eq.©. (c),(d) show the same but for k = 0.85, A = 0.47. 
All panels are for Vo = -0.5 and V = 0.019. 



wave ([H]) but now in the potential 

Vi{x) = Vocn 2 (a;,fc) - Vcn 4 (x,k), 



(9) 



and require the wave amplitude to be Ao = \/Vq + k' 2 
and frequency to be uj = 1/2 — k 2 . This immediately 
leads to the dissipative part of the potential in the form 



Wi(x) = - V 8ysn(x)cn(a;)dn(x). 

Now the hydrodynamic velocity is given by v(x) = 
■^ a -cn 2 (a;, k) and the phase is : 

6 = (W /2)E[am(x,k),k], 

where E(x, k) is the incomplete elliptic integral of the 
second kind. The obtained solutions can be either stable 
or unstable, the potential and field distributions for a 
stable solutions are shown in FigJH 

Let us now turn to the case of defocusing medium: 
a = 1 and concentrate on the conservative potential in 
the form Vi — Vbsn 2 (x, k). A stable wave can be chosen 
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FIG. 5: (color online) (a) The distributions of the real (solid 
black) and imaginary (dashed red) parts of the linear potential 
given by Eq. lfTTj) for k = 1, Ao = 0.98, (b) the distribution of 
the amplitude \ip | (solid black, left vertical axis) and the phase 
gradient v (dashed red, right vertical axis) of the field tp given 
by gO]). (c),(d) show the same but for k = 0.8, A = 1.43. 
All panels are for V ,i = 2 and W Q i = 0.4042. 




FIG. 6: (color online) (a) The distributions of the real (solid 
black) and imaginary (dashed red) parts of the nonlinear po- 
tential given by Eq. (|12|) . (|13[) . (b) the distributions of the am- 
plitude \ip\ (solid black, left vertical axis) and the phase gra- 
dient (dashed red, right vertical axis) of the field ip. (c) The 
instability developing on the solution (b) perturbed with small 
noise. All panels are for k = 0.689. 



in the form 



A = v4 dn(x,k) 



(10) 



and is induced by the spatial distributions of the dissipa- 
tion and losses 



Wi = — Wosn(a;, fc)cn(a;, k), 



9k 2 . 
-(Vo 



k 2 Al) 



(11) 



One can see that is solution is of Type 1. The frequency 
of the solutions and the hydrodynamic velocity are given 
by u = V /k 2 - 1 + k 2 /2 and v(x) = fj^dnfc). The 
phase then is readily obtained from the velocity as 9(x) = 
f_ v(y)dy. FigfS] illustrates the potential and the field 
distribution for a stable solution of such a kind. 



IV. CASE OF NONLINEAR COMPLEX 
POTENTIALS 

Let us now turn to the case where the linear potential 
is absent, i.e. Ui{x) = 0, and consider spatial dependent 
nonlincarity. Respectively it is considered that V n i{x) 
is given and the problem to design the dissipative term 
inducing the given field pattern is posed. Notice that 
in general the conservative nonlinear part cannot change 
sign, i.e. the medium cannot be focusing in some regions 
of the domain and defocusing in others. We consider 
the case of a defocusing nonlincarity a = 1 and look 
for a solution of the type (fTU|) with Aq = 1 (a solution 
belonging to Type 1) embedded in the potential 



1 



(x,k), 



(12) 



then the hydrodynamic velocity is given by v(x) = 
\f2k (l + sn 2 (ir, fe)) and imposing the frequency w = 
3/c 2 /2 + 1 + 1/k 2 nonlinear dissipative part is computed 
from © to be 

W nl = ^n{ X ^{x,k) {i fc2 _ 2fc2sn2(x5fc)) _ (13) 
dn (x, k) 

The profiles of the nonlinear potential, the distribution 
of the field and the development of the instability for the 
solution (fT6| are shown in Fig El 

In the case of focusing nonlincarity a = — 1 we con- 
sider the mode © (belonging to Type 3) loaded in the 
potential 



1 + k 2 

V n i(x) = 1 - 2 + V ni cn 2 (x, k) 
ZA 



and set uj = 1/2 — k 2 /2. Then, from © we get 

W n i(x) 



V8Vni sn(x)dn(a;) 



cn(x) 



(14) 



(15) 



The phase is: 9 = \/ 2V n i A E[&m(x, k), k]. 

Finally we consider a more general case when there ex- 
ists a linear real potential with complex linear and non- 
linear inhomogeneous dissipative terms. We let a = 1 
and consider the pattern (|TU1) loaded in the linear lattice 
Vi(x) = 2fc 2 sn 2 (x, k). Then requiring uj = k 2 + 1 we ob- 
tain that the pattern (|10|) of the Type 1 is induced by the 
distribution of linear and nonlinear dissipative terms as 
follows 



W n 



fc 2 sn(a;, k) 
dn 3 (x) 



Wi(x) =0^,(1), (16) 
[2asn(a;, k)dn(x, k) + kcn(x, k)) .(17) 
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FIG. 7: (color online) (a) The distributions of the real and 
imaginary parts of the linear potential given by Eq. (|16[) 
for k = 0.7 (they coincide), (b) the same for the nonlinear 
potential given by (|17[) . solid black line is for the real pard and 
dashed red line is for the imaginary part, (c) the distribution 
of the amplitude (solid black, left vertical axis) and the 
phase gradient v (dashed red, right vertical axis) of the field 
ip. (d) The relaxation of the initial field to the solution shown 
in (c). All panels are for a — 1 and k — 0.7. 



The phase is 8 = kx and the velocity field is therefore 
constant. The potentials and the field distributions for 
Eqs. (fT7i|) . ([TT)) when a = 1 are shown in Fig. [7] The solu- 
tion can be stable in its existence Vbz > k 2 with Vbz the 
amplitude of the real linear potential. The evolution is 
shown in panel (d) of Fig. [71 



V. DISSIPATIVE BREATHERS 

In the limit of infinite period (e.g. for modulus k = 1) 
we have that the potentials in Eqs. (|16|) -(fl7 |) give rise to 
a stable localized hump in the sense that any perturba- 
tion of it will relax back to the stationary state. Stable 
solutions of this type correspond to attractive centers. 
It is possible to destabilize it into limit cycles and cre- 
ate stable dissipative breathers, by analogy with the bi- 
furcation of a dissipative soliton to a pulsating soliton 
in the complex Ginzburg-Landau equation [12J (simi- 
lar studies can be done with other types of soliton so- 
lutions given above). We take this solution as initial 
condition of the complex NLS with a parabolic trap po- 
tential added to the real linear potential, e.g. we take 
Vi(x) = \V?x 2 + 2fc 2 sn 2 (x, k), keeping all other poten- 
tials the same. We use the frequency of the parabolic 
trap as a parameter to induce the transition into a limit 
cycle. For f2 small enough the extra potential will act 
as a small perturbation allowing the solution to relax to 
a nearby stationary solution of similar shape. We find, 
however, that when Q, exceeds a critical value the solution 
bifurcates from a single hump to a double humps shape 
which is stable and stationary e.g. it still corresponds to 
an attractive center. In the panel (a) of Fig. [5] we show 




FIG. 8: (a) The modulus of the wavefunction before (Q = 
0.05, continuous line) and after (Q = 0.075, dashed line) the 
bifurcation, (b) Relaxation of the position of the relative 
minimum of the density to the new stationary equilibrium, 
(c) Time evolution of the breating mode at 51 = 0.25. (d) 
The positition of the relative minimum of the density profile 
versus time during the breather oscillations. All panels are 
for the same parameters values as in Fig. 7 except for a = 2.0 
and k = 1. 



the shapes of the solutions taken before and after the 
bifurcation, while in panel (b) we depict the relaxation 
dynamics of the center of the two humps solution toward 
the new equilibrium center. As we increase the strength 
of the parabolic trap the relaxation time increases until 
it becomes a regular periodic oscillation (see panel (d)). 
The breather exists in a window in trap frequencies that, 
for the given parameters, starts at Hi « 0.2075 and ends 
at f^2 ~ 0.27174. Outside this window the solution re- 
turns to be stationary. In the panel (c) of Fig. [5] we have 
shown the time evolution of the density profile during the 
breather motion for a parameter value inside the above 
existence window. 



VI. CONCLUSIONS 

Summarizing, we have found the exact solitonic and 
nonlinear periodic solutions for the generalized NLS 
equation with complex linear and nonlinear potentials. 
We have confirmed numerically the stability of some so- 
lutions, in particular those having the PT symmetry. 
We have shown that if stable, the solutions are quite 
robust against perturbations and in this sense they are 
generic. We have investigated stable breather solutions 
which originate when a parabolic trap is switched on 
as additional real potential. Using the strength of the 
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parabolic trap as parameter we have demonstrated that 
the stationary solutions undergo a bifurcation to a limit 
cycle which correspond to a stable breather mode. The 
stability of both stationary periodic and single hump soli- 
ton solutions and breathers opens the possibility to ex- 
perimentally observe such waves both in nonlinear optics 
and in BEC with optical lattices in presence of dissipativc 
effects. 
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